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Abstract. Consider a linear Boltzmann equation posed on the Euclidian 
plane with a periodic system of circular holes and for particles moving at 
speed 1. Assuming that the holes are absorbing — i.e. that particles falling 
in a hole remain trapped there forever, we discuss the homogenization limit of 
that equation in the case where the reciprocal number of holes per unit surface 
and the length of the circumference of each hole are asymptotically equivalent 
small quantities. We show that the mass loss rate due to particles falling into 
the holes is governed by a renewal equation that involves the distribution of 
free-path lengths for the periodic Lorentz gas. In particular, it is proved that 
the total mass of the particle system at time t decays exponentially fast as 
t — > +00. This is at variance with the collisionless case discussed in [Caglioti, 
E., Golse, P., Commun. Math. Phys. 236 (2003), 199-221], where the total 
mass decays as Const. Jt as t — * +00. 



1. Introduction 

The homogenization of a transport process describing the motion of particles in 
a system of fixed obstacles — such as scatterers, or holes — leads to very different 
results according to whether the distribution of obstacles is periodic or random. 
Before describing the specific problem analyzed in the present work, we recall a few 
results recently obtained on a more complicated, and yet related problem. 

An important example of the phenomenon mentioned above is the Boltzmann- 
Grad limit of the Lorentz gas. The Lorentz gas is the dynamical system corre- 
sponding to the free motion of a single point particle in a system of fixed spherical 
obstacles, assuming that each collision of the particle with any one of the obstacles 
is purely elastic. Since the particle is not subject to any external force, we assume 
without loss of generality that its speed is 1. The Boltzmann-Grad limit is the 
scaling limit where the obstacle radius and the reciprocal number of obstacles per 
unit volume vanish in such a way that the average free path length of the particle 
between two consecutive collisions with the obstacles is of the order of unity. 

Call f(t,x,v) the particle distribution function in phase space in that scaling 
limit — in other words, the probability that the particle be located in an infinites- 
imal volume dx around the position x with direction in an infinitesimal element of 
solid angle dv around the direction v at time t > is f(t,x,v)dxdv. 

In the case of a random system of obstacles — more precisely, assuming that the 
obstacles centers are independent and distributed in the 3-dimensional Euclidian 
space under Poisson's law — Gallavotti proved in [121 [H] (see also [TS] on pp. 
48-55) that the average of / over obstacle configurations (i.e. the mathematical 
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expectation of /) is a solution of the linear Boltzmann equation 

(dt + v ■ X7 X + a) fit, x, v ) = — / fit, x, v — 2(ui ■ v)u>)ui ■ vdui . 



If, on the contrary, the obstacles are periodically distributed — specifically, if 
they are centered at the vertices of a cubic lattice — the limiting particle distri- 
bution function / cannot be the solution of any linear Boltzmann equation of the 
form 

(d t + v ■ \7 X + cr)f(t, x, v) = a \ p(v\w)f(t, x, w)dw , 

J\w\ = l 

where p is a continuous, symmetric transition probability density on the unit sphere: 
see [TB] for a complete proof of this negative result, based on earlier estimates on 
the distribution of free path lengths for the periodic Lorentz gas [H [T7] . 

The correct limiting equation for the Boltzmann-Grad limit of the periodic 
Lorentz gas was found only very recently: see [3 [52]. In the 2-dimensional case, 
the most striking feature of the theory presented in these references, is that the 
limiting equation is set on an extended phase space involving not only the particle 
position x and direction v, as in all classical kinetic models, but also the (rescaled) 
distance s to the next collision point with the obstacles and the impact parameter 
h at this next collision point. 

The particle motion is described in terms of its distribution function in this 
extended phase space, F = F(t, x, v, s, h), which is governed by an equation of the 
form 

(d t +v-V x - d s )F(t, x, v, s,h) = J P(s, h\h')F(t, x, R[ir - 2 arcsin(/i')R 0, h')dti 

where R[8] designates the rotation of an angle 9, and P(s,h\h') is a nonnegative 
integral kernel whose explicit expression is given in [5] but is of little interest for the 
present discussion. The particle distribution function in the classical phase space 
of kinetic theory is recovered in terms of F by the following formula: 

r + co i-l 

f(t,x,v)= / / F(t 7 x,v,s,h)dhds . 

However, the particle distribution function / itself does not satisfy a linear Boltz- 
mann equation in closed form. 

Loosely speaking, in the case of a periodic distribution of obstacles, the parti- 
cle "feels" the correlations between the obstacles, since its trajectory consists of 
segments of maximal length avoiding the obstacles. This explains the need for an 
extended phase space in order to describe the Boltzmann-Grad limit of the Lorentz 
gas, in the periodic case. In the random case studied by Gallavotti, the obstacles 
centers are assumed to be independent, which reduces the complexity of the limiting 
dynamics. 

In the present work, we shall study a much simpler homogenization problem, 
which can be formulated as follows: 

Problem. Consider a system of point particles whose distribution function is 
governed by a linear Boltzmann equation. The particles are assumed to move in 
a periodic system of holes. Describe the asymptotic behavior of the total mass of 
the particle system in the long time limit, assuming that the radius of the holes 
and their reciprocal number per unit volume vanish so that the average distance 
between holes is of the order of 1. 

Although the underlying dynamics in this problem is a lot simpler than that of 
the Lorentz gas, the homogenized equation is also set on an extended phase space, 
analogous to the one described above. 
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A we shall see, the mathematical derivation of the homogenized equation in 
the extended phase space for the problem above involves only very elementary 
arguments from functional analysis — at variance with the case of the Boltzmann- 
Grad limit of the Lorentz gas, which requires a fairly detailed knowledge of particle 
trajectories. 

2. The model 
We consider the monokinetic, linear Boltzmann equation 

(1) d t f s +vV x f e + a(f e -Kf e ) = 

in space dimension 2. 

The unknown function f(t, x, v) is the density at time t £ 1R + of particles with 
velocity v £ S 1 , located at x £ R 2 . For each <fi £ i 2 (§ 1 ), we denote 

K(f>(v) := — / k(v,w)(f>(w)dw, 

where dw is the uniform measure (arc length) on the unit circle S 1 . We henceforth 
assume that 

k £ L 2 (S 1 x S 1 ) . k(v, w) = k[w, v)>0 a.e. mv,w £ S 1 

(2) f 

and 27 / k(v,w)dw = 1 a.e. in v £ § . 

The case of isotropic scattering, where k is a constant, is a classical model in 
the context of Radiative Transfer. Likewise, the case of Thomson scattering in 
Radiative Transfer involves the integral kernel 

k(v,w) = ^(l + (vwf) 

— see for instance chapter I, §16 of [Sj. Finally, the collision frequency is a constant 
a > 0. 

The linear Boltzmann equation fTJ) is set on the spatial domain Z e , i.e. the space 
1R 2 with a periodic system of holes removed: 

Z e := {x £ JR 2 dist(a;,£Z 2 ) > e 2 } . 

We assume an absorption boundary condition on dZ e : 

f E = for (t, x, v) £ x dZ e x S 1 , whenever v ■ n x > , 

where n x denotes the inward unit normal vector to Z £ at the point x £ dZ £ . This 
condition means that a particle falling into any one of the holes remain there forever. 

The same problem could of course be considered in any space dimension. Notice 
however that, in space dimension N > 2, the appropriate scaling, analogous to 
the one considered here, would be to consider holes of radius e"'^" 1 ' centered 
at the points of the cubic lattice e7L N — see for instance [SJ [T7]. Most of the 
arguments considered in the present paper can be adapted without change to the 
higher dimensional case, except that the expression of one particular coefficient 
appearing in the homogenized equation is not yet known explicitly at the time of 
this writing. 

The most natural question related to the dynamics of the system above is the 
asymptotic behavior of the total mass of the particle system in the small obstacle 
radius e< 1 and long time limit. 

The last two authors have considered in [7] the non-collisional case (a = 0) and 
proved that, in the limit as e — > + , the solution f e converges in L°°(1R + x R 2 x S 1 ) 
weak-* to a solution / of the following non-autonomous equation: 

(3) d t f + vV x f = ^f, 
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where p is a positive decreasing function defined below. In that case, the total mass 
of the particle system decays like Const. /t sst^ +00. 

Observe that, starting from the free transport equation, we obtain a non-autono- 
mous (in time) equation in the small e limit. In particular, the solution of equation 
Q cannot be given by a semigroup in a function space such as L P (IR, 2 x S*). As 
we shall see, the homogenization of the linear Boltzmann equation in the collisional 
case (<r > 0) leads to an even more spectacular change of structure in the equivalent 
equation obtained in the limit. 

The work of the last two authors [7] relies upon an explicit computation of the 
solution of the free transport equation, where the effect of the system of holes is 
handled with continued fraction techniques. In the present paper, we investigate the 
analogous homogenization problem in the collisional case (a > 0). As we shall see, 
there is no explicit representation formula for the solution of the linear Boltzmann 
equation, other than the one based on the transport process, a particular stochastic 
process, defined for example in |23j . 

This representation formula was used in a previous work of the first author [5] , 
who established a uniform in e upper bound for the total mass of the particle system 
by a quantity of the form Const. e _a "* for some a a > 0. This exponential decay 
is quite remarkable: indeed, there is a "phase transition" between the collisionless 
case in which the total mass decays algebraically as t — » +00, and the collisional 
case in which the total mass decays at least exponentially fast in that same limit. 

In the present paper, we further investigate this phenomenon and show that the 
exponential decay estimate found in [2] is sharp, by giving an asymptotic equivalent 
of the total mass of the particle system in the small e limit as t — > +00. 

Instead of the semi-explicit representation formula by the transport process, our 
argument is based on the very special structure of the homogenized problem. The 
key observation in the present work is that this homogenized problem involves a 
renewal equation, for which exponential decay is a classical result that can be found 
in classical monographs such as |12j . 

3. The main results 

First, we recall the definition of the free path length in the direction v for a 
particle starting from x in Z e : 



The distribution of free path length has been studied in 0[T71[71|3]. In particular, 
it is proved that, for each arc IcS 1 and each t > 0, one has 



as e — > + , where |7| denotes the length of I and the measure considered in the 
statement above is the uniform measure on [0, l] 2 x S 1 . 

The following estimate for p can be found in [§]: there exist C, C > such that, 
for all t > 1: 



(4) 



T e (x,v) :=mf{t>0\x-tve dZ e } . 



(5) 



meas({0, v) e {Z e n [0, l] 2 ) x I \ er E (x, v) > t}) -> p(t)\I 




a 



uniformly as e — > + , so that 




In [3] F. Boca and A. Zaharescu have obtained an explicit formula for p: 



(8) 
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FIGURE 1. The graphs of T (left) and ofp (right) 



where the function T is expressed as follows: 
(9) 

1 



T(t) = 



24 



+ 2(l-i) 2 ln(l-i)-I(l-I) 2 ln|l-i| 



if * G (0, |], 
if t e (§,+oo) 



This is precisely at this point that the case of space dimension 2 differs from the 
higher dimensional case. Indeed, in space dimension higher than 2, the existence 
of the limit (J5]) has been proved in [21] , while the uniform estimate analogous to 
^ is to be found in [T7]. However, no explicit formula analogous to ([5]) is known 
in that case, at least at the time of this writing. We have chosen to treat in the 
present paper only the case of the square lattice in space dimension 2 as it is the 
only case where the limit ([SE]) is known completely. 

Throughout this paper, we assume that the initial data of (S £ ) satisfies the 
assumption 



(10) f " > on R 2 x S 1 and / / f n (x, v)dxdv + sup f ln (x, v) < +oo. 

JJ]R 2 xS 1 (x,u)eR 2 xS 1 

For each < £ -C 1, let f e be the (mild) solution of the initial boundary value 
problem 

' 9tf e + v ■ V,/ e + tr(f e - Kf e ) = 0, (x, v) G Z e x S\ t > 0, 

(E e )l / £ = 0, ifvn x >0, (x,v) G dZ £ x S\ 

t f e (0,x,v) = f m (x,v), (x,v)GZ s xS 1 . 

The classical theory of the linear Boltzmann equation guarantees the existence and 
uniqueness of a mild solution f s of the problem (3 E ) satisfying 



(11) 



< / £ (t, x, v) < sup f m (x, v) a.e. on H+ x Z £ x S 1 , 

(£c,«)SR 2 xS 1 



f e (t,x,v)dxdv < // f ln (x,v)dxdv . 
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Consider next F :— F(t, s, x, v) the solution of the Cauchy problem 

' d t F + v ■ V X F + d s F = -aF + j(t A s)F, t,s>0, (x, v) G R 2 x S 1 , 

p+oo 

(S) < F(t,0,x,v) = a KF(t,s,x,v)ds, t > 0, (x, v) e R 2 x S 1 , 
Jo 

F(0,s,x,v) = <7e- as f in (x,v), s > 0, (x,v) G 1R 2 x S 1 , 

with the notation t A s := min(i, s). Notice that F is a density defined on the 
extended phase space: 

{(s,x,v)\s >0,i£ R 2 ,v G S 1 } 

involving the extra variable s, whose interpretation is given below. 

Henceforth, we shall frequently need to extend functions defined a.e. on Z e by 
inside the holes (that is, in the complement of Z e ). We therefore introduce the 
following piece of notation. 

Definition : For each function tp = tp(x) defined a.e. on Z e , we denote 



tp(x) if x G Z £ 
if x $ ~Z~ e 



We use the same notation {/ £ } or {F e } to designate the same extension by inside 
the holes for functions defined on cartesian products involving Z £ as one of their 
factors, such as R+ x Z e x S 1 in the case of f £ , and R + x R + x Z e x S 1 in the case 
of F e . 

Our first main main result is 
Theorem 1. Under the assumptions above, 

r+oo 

{fe} ^ / Fds 
JO 

in L°°(R + x R 2 x S 1 ) weak-* as e — > + , where F is the unique (mild) solution of 

Notice that the limit of the (extended) distribution function of the particle 
system is indeed defined in terms of the solution F of the homogenized integro- 
differential equation (E). However, it does not seem that the limit of itself 
satisfies any natural equation. 

Next we discuss the asymptotic decay as t +oo of the total mass of the 
particle system in the homogenization limit e«l. Obviously, the particle system 
loses mass due to particles falling into the holes. 

In order to do so, we introduce the quantity: 

m(t, s) := // F(t,s,x,v)dxdv . 

A key observation in our work is that m is the solution of a renewal type PDE, as 
explained in the next proposition. 



Proposition 1. Denote 



B(t,s) = a--(tAs), 
P 



and assume that f m satisfies the condition \10\) . 
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Then the renewal PDE 

' dtfJ,{t,s)+d a fi(t,8) + B(t,s)fj,(t,s) = 0, t,s>0, 

r + co 

n(t,0)=cr / n(t,s)ds, t>0, 
Jo 

At(0,s) = ae~ as , s > 0, 

has a unique mild solution /i £ L°°([0, T]; i 1 (R + )) for all T > 0. 



Moreover, one has 

ml £_ jj 1 = — 

2vr 



IR 2 xS 1 



a.e. m (i, s) £ R+ x 1R + . 



Renewal equations are frequently met in many different contexts. For instance 
they are used as a mathematical model in biology to study the dynamics of struc- 
tured populations. The interested reader can consult |20j or |24j for more informa- 
tion on this subject. 

Consider next the quantity: 

/•+QC i- /• r+oo 

(12) M(t):=± / // F(t,s,x,v)dxdvds = / m(t,s)ds. 
* Jo iii^s 1 Jo 

As explained in the theorem below, M(t) is the total mass at time t of the particle 
system in the limit as e — > + ; besides, the asymptotic behavior of M as t — > +oo 
is a consequence of the renewal PDE satisfied by the function (i, s) <—»■ m(t, s). 

Theorem 2. Under the same assumptions as in theorem^ 
(1) the total mass 



j_ 

271 



f £ (t,x,v)dxdv -> M(t) 



in Lj oc (M,+ ) as e — » + , and a.e. in t > after extracting a subsequence of 
e -> 0+; 

(2) i/ie limiting total mass is given by the representation formula 

M (t) = 2b II f m (x,v)dxdv V K* n (t), t > 

«(t) := cre _CTt p(i)lt>o j k*" := « * . • • * « 

V v ' 

n factors 

and * denoting as usual the convolution product on the real line; 

(3) for each a > 0, there exists £ CT £ (— a, 0) such that 

M(t) ~ C^e 5 "' as i — > +00 

with 

If f n (x,v)dxdv 

1 JJlRixS 1 



2lT(T 

tp^e-^+^dt 



(4) finally, the exponential mass loss rate £ CT satisfies 
£ CT ~ — cr as a — > + , and £ CT — ► —2 as <7 ■ 
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Statement (1) above means that M is the limiting mass of the particle system 
at time t as e — > + . Statement (3) gives a precise asymptotic equivalent of M{t) 
as t — > +oo. 

As recalled in the previous section, if a = in the linear Boltzmann equation 
(S e ), the total mass of the particle system in the vanishing e limit is asymptotically 
equivalent to 



as t — > +oo. The reason for this slow, algebraic decay is the existence of channels 
- infinite open strips included in the spatial domain Z e , i.e. avoiding all the 
holes. Particles located in one such channel and moving in a direction close to 
the channel's direction will not fall into a hole before exiting the channel, and this 
can take an arbitrarily long time as the particles' direction approaches that of the 
channel. This construction based on channels leads to a sufficiently large fraction 
of the single-particle phase space and accounts for the algebraic lower bound in 
(J7|) . The asymptotic equivalent mentioned above in the collisionless case a = is a 
consequence of a more refined analysis based on continued fractions given in [7]. 

When a > 0, particles whose distribution function solves the linear Boltzmann 
equation in (H E ) travel on trajectories whose direction is discontinuous in time — 
more specifically, time discontinuities are distributed under an exponential law of 
parameter a. Obviously, this circumstance destroys the channel structure that 
is responsible of the algebraic decay of the total mass of the particle system in 
the collisionless case, so that one expects that the total mass decay is faster than 
algebraic as t — > +oo. That this decay is indeed exponential whenever a > is by 
no means obvious: see the argument in [5J, leading to an upper bound for the total 
mass. Statement (3) above leads to an asymptotic equivalent of the total mass, 
thereby refining the conclusions of [5J ■ 

In section 4, we give the proof of theorem [TJ the evolution of the total mass in 
the vanishing e limit (governing equation and asymptotic behavior as t — ► +oo) is 
discussed in section 5. 



Our argument for the proof of Theorem [JJ is split into several steps. 

4.1. A new formulation of the transport equation. Perhaps the most surpris- 
ing feature in Theorem [1] is the introduction of the extended phase space involving 
the additional variable s. 

As a matter of fact, this additional variable s can be used already at the level 
of the original linear Boltzmann equation — i.e. in the formulation of the problem 



Let us indeed return to the initial boundary value problem (S e ) for the linear 
Boltzmann equation. 

As recalled above, the last two authors have obtained the homogenized equation 
corresponding to (S E ) in the noncollisional case (a = 0) by explicitly computing the 
solution of the linear Boltzmann equation for each < £ < 1. In the collisionnal 
case (cr > 0), as recalled above, there is no such explicit formula giving the solution 
of the linear Boltzmann equation — except the semi-explicit formula involving the 
transport process defined in [23) . 

However, not all the information in that semi-explicit formula is needed for the 
proof of Theorem [TJ The additional variable s is precisely the exact amount of 




TP 



H 



4. The homogenized kinetic equation 
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information contained in that semi-explicit formula needed in the description of 
the homogenized process in the limit as e — > + . 

Consider therefore the initial boundary value problem 

d t F £ + v ■ V X F £ + d s F £ + oF £ = 0, t,s>0, (x, v) e Z e x §\ 



F e (t, s, x, v) — 0, if v ■ n x > 0, t, s > 0, (x, v) 6 (dZ £ x 8 1 ), 

/>oo 

F e (t,0,x,v) = a KF £ (t,s,x,v)ds, t > 0, (x, v) e Z £ x S 1 , 
Jo 

s > Q,(x,v) £ Z £ x S 1 , 



[ F e (0, =ae— / m (x, W ), 

with unknown i 7 ^ := F £ {t, s, X, v). 

The relation between these two initial boundary value problems, (H e ) and (S £ ), 
is explained by the following proposition. 



Proposition 2. Assume that f m satisfies the assumption U0\) . Then 
a) for each e > 0, the problem (S E ) has a unique mild solution such that 



(t,x,v) 



\F e (t,s,x,v)\ds belongs to L°°([0,T] x Z e x S 1 ) 



for each T > 0; 
&J moreover 

< F e (t,s,a:,i;) < ||/ iri ||ioo (IR2 xS i )C re" 
a.e. int,s>0, x £ Z £ and v £ § , and 

F e (t,s,x,v)ds = f £ (t,x,v), 



for a.e. t>0,x£Z £ ,v£§ 1 , where f £ is the solution of (3 £ ). 

Proof. Applying the method of characteristics, we see that, should a mild solution 
F £ of the problem (S £ ) exist, it must satisfy 



(13) 
with 



F £ (t,s,x,v) = Fi yS (t,s,x,v) + F2, e (t,s,x,v), 



Fi tB (t,s,x,v) = l s<eTe( |^)l s < t e as F £ (t - s,0,x - vs,v) 

(14) _ r+oo 

= ^ s <er € (2i,v)^s<tO-e as I KF £ (t — s,t,x — sv, v)dr 
Jo 



and 



(15) 



F 2 , £ (t,s,x,v) = t t <er e (^,v)h< s e CT *F E (0,s - t,x- vt,v) 



= h< er ^,v)U<s(re-' , 'f m (x-tu,v) 

a.e. in (i, s, x, v) £ R + x 1R + x R 2 x S 1 . 

First, define Xt to be, for each T > 0, the set of measurable functions G defined 
on E + x Ki x Z £ x S 1 such that 



(t,x,v) 



r+oo 

/ \G(t,s,x,v)\ds belongs to L°°([0,T] x Z £ x S 1 ) 
Jo 



which is a Banach space for the norm 



l|G|| 



\G(;S,;-)\ds 







L°°([0,T]xZ e xS!) 
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Next, for each G € Xt, wc define 



r+ao 

TG(t,s,x,v) := t s <er e ^, v )^s<t<re-" s / KG(t 

Jo 



s,t,x — sv, v)dr . 



Obviously 



L°°(Z e xS 1 ) 



r-l-oo 

/ \T n G(t,s,-,-)\ds 
Jo 

pt r + co 

<a / \T n - 1 G(t 1) T,',')\dT 
Jo Jo 



'0 Jo 



L°°(Z e xS 1 ) 



+ 00 



\G(t n ,s, -,-)\ds 



dt n . . . dt\ 



L°°(Z e xS!) 



(aT) n 



so that 



Now Fi t e = 1~F S , so that (|T3|) can be recast as 

F £ = F 2i£ + TF £ . 

This integral equation has a solution F £ S Xt for each T > 0, given by the series 

F £ = T n F2,e 



which is normally convergent in the Banach space Xt since 

E \\ Tnp ^\u T < E ^f\\^e\\ XT < +™ ■ 

n>0 n>0 

Assuming that the integral equation above has another solution F £ 6 Xt would 
imply that 



F s - F' e = T(F £ - F e ') = . . . = r T n {F £ - F' e ) , 



so that 



\T n {F e -F' e )\\ XT < 



(aT) n 



F K -Ft\ 



At 







as n — > +00: hence F' e — F e . Thus we have proved statement a). 

As for statement b), observe that TG > a.e. on R + x M + x Z £ x S 1 if G > 
a.e. on K+xK+x^xS 1 . Hence, if f in e L°°(R 2 x S 1 ) satisfies f in > a.e. 
on R 2 x S 1 , one has F 2 , e > a.e. on 1R+ x R + x Z £ x §\ so that T n F 2 , £ > 
a.e. on 1R + x 1R + x Z £ x S 1 and the series defining F £ is a.e. nonnegative on 
R+ x R+ x Z £ x S 1 . 

Next, integrating both sides of (fTT?|) with respect to s, and setting 

r+00 

g E (t,x,v):= F e (t,s,x,v)ds , 
Jo 
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we arrive at 

r+oc r + co 

g e (t,x,v) = / F 2 . e (t, s,x, v)ds + / Fi i£ (t,s,x,v)ds 
Jo Jo 

r+oc 

= l t <er E (^,v)f n (x - tv,v) / l t<s <Je~ as ds 

Jo 

r+oo / I-+00 \ 

+ / l s<eTe (|^)l s <tcre _<TS f / KF e (t~s,r, x-sv,v)dT\ ds 

= l t <er e (f,v)f m (x-tv,v)e- at 

+ / e~ as l s<£Tr (2L tV )oKg e (t - s,x - sv,v)ds 
Jo 

in which we recognize the Duhamel formula giving the unique mild solution f e of 
(S E ). Hence 

/•+0O 

f e (t,x,v) = / F e (t, s, x, v )ds a.e. in (t, x, v) S 1R + x Z £ x S 1 . 



Finally, since (^, e ) satisfies the maximum principle, one has 

f £ (t,x,v) < ||/ m ||L-(R2 X si) a.e. in (t,x,v) E R + x Z e x S 1 . 
Going back to (fll?)) . we recast it in the form 

F e (t,s,x,v) = t s<eT ^^)l s<t (Te^' ys K f £ (t - s,x - sv,v) 
+ lt<eT s (f,v)h<sCre- as f m (x-tv,v) 

< ls<£T e (f ,t))ls<tO'e _0 ' S ||/ m || L = ( R 2 x gi : | 
+ l*<ETe(f ,D)lt<S Cre ~ <TS |l/ m |U=(R 2 xS 1 ) 

Sere ||/ ||l=(r2 x si) 
a.e. in (t, s, x, v) S R + x R + x Z e x S , which concludes the proof. □ 
Observe that if 

F e (0,s,x,v)=ae-° s f in (x,v) 

is replaced with 

F e {Q, S ,x,v)=R(s)f in (x,v) 

where II is any probability density on R + vanishing at oo, the conclusion of the 
lemma above remains valid. In other words, the dependence of the solution F s of 
the problem (S) upon the choice of the initial probability density II disappears after 
integration in s, so that the particle distribution function f e is indeed independent 
of the choice of II. 

The extra variable s in the extended phase space has the following interpretation. 
Recall that the solution f £ of the linear Boltzmann equation can be expressed in 
terms of the transport process (see 23 ), a stochastic process involving a jump 
process in the v variable, perturbed by a drift in the x variable. The variable s is 
the "age" of the current velocity v in that process, i.e. the time since the last jump 
in the v variable. 

The choice II(s) = ae" as corresponds with the situation where the gas molecules 
have been evolving under the linear Boltzmann equation for t < and the holes 
are suddenly opened at t = 0. 

Before giving the proof of Theorem [I] we need to establish a few technical lem- 
mas. 
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4.2. The distribution of free path lengths. A straightforward consequence of 
the limit in §E§ is the following lemma, which accounts eventually for the coefficient 
p(t A s)/p(t A s) in the limiting equation (E). 

Lemma 1. Let t e be the free path length defined in Q. Then for each t > 
in L°°(R 2 x S 1 ) weak-* as e -> 0+. 

(See the definition before Theorem [T] for the notation {lt<eT t (-,v)}-) 
Proof. Since the linear span of functions x = x( x j v ) °f the form 

(/.(Xju) = x (x)ti(v), 0eC^°(R 2 ) and / an arc of S 1 

is dense in L 1 (R 2 x S 1 ), and the family l ETe (* lt! )>t is bounded in L°°(R 2 x S 1 ), it 
is enough to prove that 



Write 



4>{x, v)l STe (£ tV -) >t dxdv — ► p(t) // v)dxdv as e — > . 



(x,v)l eTE (2L tV)>t dxdv = J x{x) (^J leT e (f ,v)>tdvj dx 



>,v)>tdv . 



with 

r e (y) := ^ l £ r £ (a,t 
Obviously T £ is 1-periodic in y\ and j/2 and satisfies < T e < |J|. Hence 
l W) >eTe(y) = X! r E (A;)e 2iwfe ^ 

in L 2 (R 2 /Z 2 ) with 

f e (fc) := / T e (z)- 2lTrfe - z dz 
Jm^Qm |,|* a l)<i/2 
|*|>e 

for each k &TL 2 . 

Then, by Parseval's identity, 

^(x)T £ (-) dx = f X (x) ( ]T T e (fc)e 2 ^ ) dx 



2 e 



= x(0)f e (0) + 53 f e (fc)x(-2^fc/e) 
fcez 2 \(o,o) 

with 

Jr. 2 

Applying again Parseval's identity, 

E i^( fc )i 2 = / (l , I5<1/ J r '(v)i 3 */<|J| 

/ max(| yi I ,|j/2 1X1/2 
fceZ 2 |B|>£ 

while 

lx(OI<^il|v 2 xll^, 
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so that 



2 



|*(-27rfc/e)| < j-^W^X 



4tt 2 |£| 21 

Hence, by the Cauchy-Schwarz inequality, 

2 



fcez 2 \(o : o) 



fcez 2 \(o,o) fcez 2 \(o,o) 1 1 



and therefore 

X (:r)T £ (£) = x(O)XUO) + 0(e 2 ) 



as e — * + . 
By © 



so that 



t e (0)= f T s (y)dy ^ p(t)\I\ as £ ^0+. 

IhI>c 



X(0)T E (0) ->p(i)|/| / x(x)dx=p(t) cb{x,v)dxdv 
Jm. 2 JJ-RixS 1 

as e — > + , and hence 

x(x)T e (-) da; = / / v)dxdv + o(l) + 0(e 2 ) 

which entails the announced result. □ 

4.3. Extending f e by in the holes. We begin with the equation satisfied by 
the (extension by inside the holes of the) distribution function {/ e }. 

Lemma 2. For each e > 0, the function {f e } satisfies 

(d t +v V x ) {f e } + a({f E } - K {/J) = (v ■ n x )f e \ dZeXS1 5 dZe 

in D'QRi xR 2 xS 1 ), where Sq Zs is the surface measure concentrated on the boundary 
of Z £ . and n x is the unit normal vector at x £ dZ e pointing towards the interior of 
Z e . 

Proof. One has 

dt {fe} = {dtfe} 

and 

Va; {f e } = {VJ e } + f e \dZ e xS^ 5 dZs n x 

in T>'(JR* + x M 2 x S 1 ). Hence 

= {d t f s + v-V x f s + a(f e -Kf s )} 

= d t {f £ } + v-V x {f £ } + (v ■ n x )f e \ gZeySl 6 dZc + a({f £ } - K {/ e }) 

in V'(R* + x K 2 x S 1 ). □ 

A straightforward consequence of the scaling considered here is that the family 
of Radon measures 

{v-n x )f E \ d ^ xS1 S dZe 
is controlled uniformly as e — > + , in the following manner. 
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Lemma 3. For each R > 0, the family of Radon measures 



(v ■ n x)fe\ aZeyM i&dz E \ [ _ RB ^ yM i 
is bounded irfl M([-R, R} 2 X S 1 ). 



Proof. The total mass of the measure 



is less than or equal to 

27r||/ e ||L~(]R + xz e xsi)||'5az e \[-r,r] 2 \\m([-rm] 2 ) 
which is itself less than or equal to 




as e — > + , whence the announced result. 



□ 



4.4. The velocity averaging lemmas. As is the case of all homogenization re- 
sults, the proof of Theorem [T] is based on the strong L\ convergence of certain 
quantities defined in terms of F £ . In the case of kinetic models, strong Lj com- 
pactness is usually obtained by velocity averaging — see for instance [TJ QUI [18] 
for the first results in this direction. Below, we recall a classical result in velocity 
averaging that is a special case of theorem 1.8 in [5], 



Proposition 3. Let p > 1 and assume that f e = f e (t,x,v) is a bounded family in 
LLJJ&t xB^ xSf 1 ) such that 



for each T > and R > 0. Then, for each i/j e C(§ d_:L x S d_1 ), the family p^[f e ], 
defined by 



is relatively compact in i ; 1 oc (R t F x x Sjj 1 ). 

A straightforward consequence of Proposition [3] is the following compactness 
result in Lj oc strong, which is the key argument in the proof of Theorem [T] 

Lemma 4. Let f e = f e (t,x,v) be the family of solutions of the initial boundary 
value problem (H e ). Then the families 



are relatively compact in L ; 1 oc (R+ x R 2 x S 1 ) strong. 



For each compact subset K of H , we denote by M(K) the space of signed Radon measures 
on K, i.e. the set of all real-valued continuous linear functionals on C(K) endowed with the 
topology of uniform convergence on K. 





K{f e } = {Kf e } 



and 
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Proof. We recall that, by the Maximum Principle for (E e ), 

\fe(t,X,v)\ < ||/ m || L ^ (R2xS1) 

a.e. in t > 0, x G Z e and v € S 1 , so that 

(16) SUp || {/ e } ||l~(R+x]R 2 xS1) < ||/*" (R 2 xS 1 ) ■ 

e 

By Lemma [51 satisfies the equation 

dt {/J + v ■ V x {/J = a(K {f s } - {/,}) - S BZ . (v.n x )f £ | eZeXSl 

in 2?'(R5_ x R 2 x S 1 ). Because of (fT6]) and the fact that the scattering kernel k is 
a.e. nonnegative (see ©), one has 

\W(K{fe} ~ {A})IU-(R+x]R 2 xSi) < cr(l + ||ifl|| L =c (S i))|| {/J || L » (]R+ xR 2 xS i ) 

= 2(J II {/e} |U=°(R + xR 3 xS 1 ) 

since XI = 1 (see again ([2]).) Besides the family of Radon measures 

Ms = fs |az e x§i (« • n x )8 d z e 

satisfies 

sup / |/x e | < +oo 

e J[0,T]xS(0,fl)xS 1 

for each T > and i? > according to lemma [3J 

Applying the Velocity Averaging result recalled above implies that the family 

g e dv 

s 1 

is relatively compact in L] oc (H + x R 2 xS 1 ). 

By density of C(S 1 x S 1 ) in i 2 (§ 1 xS 1 ). replacing the integral kernel k with a 
continuous approximant and applying the Velocity Averaging Proposition [3] in the 
same way as above, we conclude that the family Kg £ is also relatively compact in 
Ll oc (R+ x R 2 x S 1 ). □ 

4.5. Uniqueness for the homogenized equation. Consider the Cauchy prob- 
lem with unknown G = G(i, s, x, v) 

(fit + v-V x + d s )G = -oG + *y AS ? G, f,s>0,ie R 2 ,v € S 1 , 

p(i A sj 

G(t,0,x,u) = S(t,x,v), t > 0, (ar,u) G R 2 x S 1 , 

G(0, s, ar, «;) = G™(s, a, u), s > 0, (x, v) G R 2 x S 1 . 

If, for a.e. (t, s, x, G R+ x R + x R 2 x S 1 , the function r G(i+r, s+r, x + ru, v) 
is G 1 in t > 0, then, since the function p G C 1 (R+) and p > on R + , one has 

T" + CT _ ^77 — G (* + T, * + t, s + rw, t> 
rfr p(i A s + r) / 

= e ffT p(t A s + t)— j 2 — = . 

FK ' dr \ p(tAs + r) J 

Hence 

e aT G(t +t,s + t,x+ tv, v) 

I : th 

p(t A s + t) 

is a constant. Therefore 

T(-t) if t < s, 



r(o) 



r(-s) if s < t, 
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SO that 

G(t,s,x,v) = l t<s e~ at p{t)G m (s -t,x- tv,v) + l s<t e" as p(s)S(t -s,x- sv,v) . 

Proposition 4. Assume that f™ £ L°°(R 2 x S 1 ). Then the problem (E) has a 
unique mild solution F such that 

r+oo 

(t, x,v)^ / \F(t, s, x, v)\ds belongs to L°°([0, T] x R 2 x S 1 ) 
Jo 

for each T > 0. This solution satisfies 

F(t, s, x, v) = l t<s ae- at p{t)f m (x - tv, v) 

r+oo 



r+oc 

+ l s<t o-e~ crs p(s) / KF(t - s,t,x - sv,v)dr 
Jo 



for a.e. (t, s, x, v) £ R + x R + x R 2 x S 1 . 

Besides, F > a.e. on R + x R + x R 2 x S 1 if f 171 > a.e. on R 2 x S 1 . 

Proof. That a mild solution of the problem (E), should it exist, satisfies the integral 
equation above follows from the computation presented before the proposition. 

As above, let be, for each T > 0, the set of measurable functions G defined 
a.e. on R + x R + x R 2 x S 1 and such that 



(t,x,v) 



+ 00 







which is a Banach space for the norm 

f + oo 

\y T 



\G{t, s, x, v)\ds belongs to L°°([0, T] x R 2 x S 1 ) , 
the norm 

r + co 

/ \G(;S,;-)\ds 

Jo 



Next, for each G £ yr, we define 



QG(t, s, x, v) := l s<t ae-° s p(s) / KG(t 

Jo 



s,t,x — sv, v)dr . 



Since < e~ as p(s) < 1, the integral kernel k > on S 1 x S 1 and Kl = 1 by ©, 
one has 



ft p + oo 



/ \QG(t,s,x,v)\ds <a 
Jo Jo Jo 

a.e. in (t, x, v) £ [0, T] x R 2 x S 1 , meaning that 

r + oo 

/ \Q n G(t,s,;-)\ds 

Jo L=°(R 2 xS!) 

pt r + co 

<<T / IQ^G^S,;-)^ 

Jo JO 



G(t-S,T,;-)\dT 



ds 



dti 



L°°(R 2 xS 1 ) 



< 



ft rt n -i r+oo 

o- n / ... / / \G(t n ,s,-,-)\ds 

Jo Jo JO 



dt n . . . dti 



L°'(R 2 xS 1 ) 



In particular 



\Q n G\ 



3t 



< 



(ory 



\G\ 



The integral equation in the statement of the proposition is 

F = F 2 + QF 

where 

F 2 (t,s,x,v) = l t<s ae- at p(t)f m (x-tv,v). 
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Therefore, arguing as in the proof of Proposition one obtains a mild solution of 
(£) as the sum of the series 

F = J2 Q n F 2 , 

n>0 

which is normally convergent in the Banach space for each T > 0. 
Should there exist another mild solution, say F', it would satisfy 

(F - F') = Q(F - F') = ... = Q n (F - F') 

for all ft > 0, so that 

\\F - F'\\y T = \\Q n (F - F')\\y T < { ^\\F F'\\y T -> 

nl 

as n — > +oo, which implies that F = F' a.e. on R + x R + x R 2 x S 1 . 

Finally, QF > a.e. on R+ x R + x R 2 x S 1 if F > a.e. on R+ x R + x R 2 x S 1 . 
Since F is given by the series above, one has F > a.e. on R + x R + x R 2 x S 1 
whenever f m > a.e. on R 2 x S 1 . □ 



4.6. Proof of the homogenization theorem. Start from the decomposition (fT3|) 
of F e . Passing to the limit as e — > + in the term i 7 ^ is easy. Indeed, by Lemma [T] 

(17) {W.(*,„)}^P(*) 

in L°°(R 2 x Sj) weak-* for each t > 0, as e — > + . Hence 

{^2,e}(i,S,X,7j) =l t<se - CTS / m (^-to,f){W e (f,-)} 

1 LX ' ^ lt <s e-° s f n {x - tv, v)p(t) =: F 2 (t, s, x, «) 

in L°°(R+ x R+ x R 2 x Sj) weak-* as e -> 0+. 

Next, we analyze the term Fx >e ; this is obviously more difficult as this term 
depends on the (unknown) solution F £ itself. 

We recall the uniform bound 

SU P II {fe} Ili-tlR+xR^S 1 ) < HZ' "H-t^OR^xS 1 ) 

— see Proposition [2] b) , so that, by the Banach- Alaoglu theorem 
(19) {f £ } / in L°°(R + x R 2 x S 1 ) weak-* 

for some / 6 L°°(R+ x R 2 x S 1 ), possibly after extracting a subsequence of e — » + . 
Thus, applying the strong compactness Lemma |4] shows that 

if {/ £ } -» if/ in Lj oc (R, + x R 2 x S 1 ) strong 

as e — > + . 

This and the weak-* convergence in Lemma [T] imply that 

{*i, e } =l s<t ae- ffS if {/J (t - s,x - v)l s<STe(f >o) 
l s<t <re~ c ' s Kf(t - s, x - sv, v)p(s) 
in L[ oc (R+ x R + x R 2 x S 1 ) weak as e -> 0+. Therefore 

{F e } (t, s, x, v) ^t s<t cre~ as Kf{t - s, x — sv, v)p(s) + F 2 {t, s, x, v) 
= : F(t, s, x, v) 

in Lj oc (R + x R + x R 2 x S 1 ) weak as e -> 0+. 
Fix T > 0; then, for t G [0, T], one has 

/ F e (t,s,x,v)ds = F 1 . e (t,s,x,v)ds + e- at f n (x-tv,v)l t<eTe ^ vV) 
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since F\ >e is supported in s < t < T, so that 



roo r 1 

I {F e } (t,s,x,v)ds — v / l s <tK f(t — s, x — vs,v)ae~' TS p(s)ds 
Jo Jo 

(21) +f n (x-tv,v)e-° t p(t) 

F(t, s, x, v)ds 



in Ll oc (R + x S 2 x S 1 ) weakly as e — > + . On the other hand 

{F e } (i, s, x, v)ds = {f e } (t, x, v) f(t, x, v) 



in L°°(R + xl 2 x S 1 ) weak-* as e -> 0+ — and therefore also in i ; 1 oc (R+ x R 2 x S 1 ) 
weak as e — > + . By uniqueness of the limit, we conclude that 

/>oo 

(22) /(*, x, v) = I F(t, s, x, v)ds a.e. in (t, x, v) £ R+ x R 2 x S 1 

so that F satisfies 



(i 



F(t, s, x, v) = t s<t cre (7S K I / F(t — s,u,x — sv, -)du ) (v)p(s) 
(23) \Jo J 

+ l t<s ae-™f in (x-tv,v)p(t) 

a.e. in (t,s,x,v) £ R+ x R + x R 2 x S 1 . By Proposition 2J this means that F is a 
solution of the Cauchy problem (S). 

By uniqueness of the solution of (E), we conclude that F — F, and that the 
whole family 

F £ -± F in Ll oc (R+ x R + x R 2 x S 1 ) 

weakly as e — > + . 

Finally, (19]) and {22]) imply that 

p OO 

{fe} ^.f= Fds 
Jo 

in L°°(R + x R 2 x S 1 ) weak-* as e — > + , which concludes the proof of Theorem[TJ 

□ 



5. Asymptotic behavior of the total mass in the long time limit 

The formulation of the homogenized equation (problem (E)) as an integro- 
diffcrcntial equation set on the extended phase space involving the additional vari- 
able s is of considerable importance in understanding the asymptotic behavior of 
the total mass of the particle system as the time variable t — > +oo. Indeed, this 
formulation implies that the total mass of the particle system satisfies a renewal 
equation, i.e. a class of integral equations for which a lot is known on the asymp- 
totic behavior of the solutions in the long time limit — see for instance in [T^] the 
basic results on renewal type integral equations. 



5.1. The renewal PDE governing the mass. We begin with a proof of Propo- 
sition [T] 
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Proof. That fj, is a mild solution of the renewal PDE means that, for a.e. (t, s) £ 
R+ x R+, 

/■+OO 

s) = l t<s ae-^ s -^e- at P (t) + l s<t e- <TS p(a) / p(i - s, r)dr 



./O 

= <re -< "p(t A s) (l t<s + l s<t J n(t - s, t)o!t 
For each T > 0, define 

/•+CSO 

Hli(t, s) = l s<t ae- as p{s) / n(t - s, t)oIt 

Jo 

a.e. in (t,s) £ R + x R + . Obviously, for each £ L°°([0, T]; i 1 (R + )) and a.e. 
t > 0, 



IIW,-)IUi(R + ) < / ^^ S) p(t-s)||0( S ,-)llL 1( R + )rf S 

< G \ ll<K«> , )IU 1 (R+) ds > 



so that, for each n > 0, one has 

||ft n #,-)IUHB+) ^ / / \\<P(tn,-)\\L^n+)dtn...dt 1 



JO 

< ^*)"iuii 

^r"^" Loo([o ' TI;L1(R+)) 

a.e. in t £ R + . 

Arguing as in the proof of Proposition [51 we see that the renewal PDE has a 
unique mild solution /j, £ L°°([0, T]; i 1 (lR+)) for all T > 0, which is given by the 
series 



n>0 

where 

Obviously IZcj) > a.e. on R + x R + if > a.e. on R + x R + , so that fi > 
a.e. on R + x R + . Besides, for each T > 0, 

IMIl~([o,t];L1(R+)) < ""j - H^ m Hi 1 (R+) = e<TT ' 

which implies in turn that 

< /u(t, s) < ae~ as p{t A s) (l t<s + l s<t e CTT ) < cre^e - '™ 

a.e. in (t,s) £ [0,T] x R + . 

Finally, let F be the mild solution of the problem (E) obtained in Proposition^ 
Since F > a.e. on R + x R + x R 2 x S 1 is measurable, one can apply the Fubini 
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theorem to show that 



i(t,s) : — / / F(t,s,x,v)dxdv 

= t t<s ae- at p{t)^ f [ f m (x - tv, v)dxdv 

J J^xS 1 

+ l t<s ae~' yt p(s) / 5¥ / / KF(t - s,t,x — sv,v)dxdvdr 

JO ii^xS 1 

= l t<s ae-^p(t)^ [ [ f m {y,v)dydv 

J J R 2 xS 1 

+ l t< s<Je~ at P (s) f ± [[ KF(t - s, t, y, v)dydvdr 

JO JJ^xS 1 

'p(*)a? // f in (y,v)dydv 

J J R 2 xS 1 



lt<s 



+ l t<s cre / 2¥ / / F(t~ s,T,y,w)dydwdT 

JO JJR 2 xS 1 

= l t<s(T e- fft p(t)^ / / - to, v)dzdt> 

J J R 2 xS 1 



+ li< s ae CTt p(s) / m(t-s,T)dT, 
Jo 

where the second equality follows from the substitution y = x — tv that leaves the 
Lebesgue measure invariant, while the third equality follows from the identity 



± / k{v,w)dv = 1 , 
Js 1 

which implies that 

± / KF(t - s,T,y,v)dv = ± F(t - s,T,y,w)dw . 
Js 1 Js 1 

In other words, 

mfc s) satisfies the same integral equation as ^ [[ V )dyd V . 
Now the solution f e of (S e ) satisfies 
f £ > a.e. on R + x R 2 x 8 1 and // f £ (t,y,v)dydv <// f ln (y,v)dydv, 

JJ^xS 1 JJ^xS 1 

which implies by Theorem [T] that 

fc(t,y,v)dvdy / / F(t,s,y,v)dvdyds . 

Ivl^RJs 1 Jo Jlvl^RJs 1 

Hence, by Fatou's lemma 

F(t, s,y,v)dvdyds < lim // f E (t, x,v)dxdv 

ylKRJS 1 s^O+JJ^xS 1 

< [[ f m (y,v)dydv, 

a.e. in t > 0. 

Letting R — > +og in the inequality above, we see that to e L°°(R + ;L 1 (R+)) 
and we have proved that the difference 

A(i, a) = TO(t, fl ) - f /" r[yi v)dydv 

ZlT J v/R 2 xS 1 
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satisfies 

A G L 00 (R + ;L 1 (R+)) and A = TZA . 

By the same uniqueness argument as in the proof of Proposition [?J we conclude 
that A = a.e. on R + x R + . □ 

5.2. The total mass in the vanishing e limit. By Theorem [TJ the solution f E 
of (S e ) satisfies 

r + oo 

{fe} / Fds in L°°(R + x R 2 x S 1 ) weak-*; 
Jo 

therefore, checking that 

p p P + OO p p 

/ / {f E }dxdv / / / Fdxdvds =: 2nM(t) 

reduces to proving that there is no mass loss at infinity in the x variable. 
Lemma 5. Under the same assumptions as in Theorem]^ 

± [[ f E (t,x,v)dxdv=± [I {f E }(t,x,v)dxdv -» M(t) 
JJZeXS 1 JJtr, 2 xs 1 

strongly in L ; 1 oc (R + ) as e — > + . 

Proof. Going back to the proof of Proposition [2] (whose notations are kept in the 
present discussion), we have seen that 

F e = rnF 2,e on R+ x R + x Z e x S 1 , 
n>0 

with the notation 

F 2 ,s{t, s,x,v) = l t<ETe (*L tV )t t < s o-e~ as f ln (x - tv,v) . 
Since T$ > a.e. whenever $ > a.e., the formula above implies that 
F E < G := Y r " G 2 a.e. in (t, s, x, v) G R+ x R + x Z E x S 1 , 

Tl>0 

where 

G 2 (t,s,x,v) := l t<s ae-° s f in {x-tv,v). 
Thus, G satisfies the integral equation 

G = G 2 + TG 
meaning that G is the mild solution of 

' {d t +v-y x + d s )G = -aG , t, s > , x G R 2 , \v\ = 1 , 

r + oo 

G(t, 0, x, v) = a I KG{t,s,x,v)ds , t > 0, x G R 2 , \v\ = 1 , 
Jo 

G(0, s, sc, v) = f m (x, v)ae- as , s > , i G R 2 , \v\ = 1 , 

Reasoning as in Proposition [2] shows that 

r + OO 

g(t, x, v) := / G(t 7 s,x,v)ds 
Jo 

is the solution of the linear Boltzmann equation 

(9 4 + v ■ V x )g + a(g-Kg) = 0, t > , x G R 2 , \v\ = 1 , 

ff (0,x, W ) = / m (x, W ), i£E 2 , |«| = 1. 
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In view of the assumption (TTU)) bearing on f m , we know that 
G > a.e. on R + x R+ x R 2 x S 1 

and 

/ // G(t,s,x,v)dxdvds = // g(t 1 x 1 v)dxdv 
Jo JJ^xS 1 JJrzxS 1 

/ m (x, v)dxdv 

for each t > 0. 

Summarizing, we have 

< {F £ } < G 

and 

/ / / G(t, s,x,v)dsdxdv = f m (x,v)dxdv < +oo . 

J J Jr+x^xS 1 JJ^xS 1 

Then we conclude as follows: for each R > 0, one has 

f E (t,x,v)dxdv — / // F(t, s,x,v)dxdvds 

Z^xS 1 JO JJ^xS 1 

{F e }(t, s, x, v)dvdxds 

\x\yRJs 1 

/ / / ({F e } — F) (t, s, x, v)dvdxds 
Jo JmkrJs, 1 

{F}(t,s,x,v)dvdxds = In, e (t) + II R>e (t) + HI R (t) . 



x\>R JS 

First, for a.e. t > 0, the term lR. £ (t) — > as R — > +oo uniformly in e > since 
< {F s } < G and G 6 L°°(R + ; x R 2 x S 1 )). 

Next, the term II R ^(t) — > strongly in L i 1 oc (R + ) as e — > + for each R > by 
Lemma SJ 

Finally, since {F e } — F in i ; 1 oc (R + x R + x R 2 x S 1 ) weak as e — > + , one 
has < {F} < G, so that F € J L°°(R + ; i 1 (R+ x R 2 x S 1 )). Hence the term 
III R (t) -> as i? -> +oo for a.e. i > 0. 

Thus we have proved that 

f e (t,x,v)dxdv — > / // F(t, s, x,v)dxdvds 



in L ; 1 oc (R + ) and therefore for a.e. t > 0, possibly after extraction of a subsequence 
ofe^0+. □ 

5.3. An integral equation for M. Given a function ip defined (a.e.) on the 
half-line R+, we abuse the notation - 01r + to designate its extension by on M,*_. 
Henceforth we also denote 

:= p(t)ae- at l t > - 

Lemma 6. The function M defined in hl l 2fy satisfies the integral equation 

M(t) = k * (Ml n+ )(t) + 5^K(t) / / f in (x,v)dxdv, t > 
where * denotes the convolution on the real line. 
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Proof. We apply the same method as for deriving the explicit representation formula 
for F starting from the equation in Corollary [H in order to find an exact formula 
for m. Indeed, by the method of characteristics, 

m(t,s) = l s< tp{s)e- as m{t- s,0) + l t<sP (t)e-' Jt m(0,s-t) 

poo 

= t s< tp{s)ae~ as / m(t — s,u)du 
Jo 

+ l t<s p(t)o-e-° s ± ff f m (x,v)dxdv. 



2w , , 

J JH 2 xS L 

The function m satisfies therefore 

m(t, s) = l s<t p(s)ae~ (TS M(t - s) 

<24) +l t<s p(t)ae-° s ±- ff f m (x,v)dxdv. 

We next integrate both sides of ([2"4")) in s E R + . By the definition (JT3J) of M, we 
obtain 

M(t)= f ap{s)e-" JS M{t- s)ds + p(t)e~ at ± I f f ln {x,v)dxdv 

Jo ^JJ^xS 1 

a.e. in t > 0, which is precisely the desired integral equation for M: 
(25) M(t)= f K(s)M(t- s)ds+ 2^n(t) f f f n (x,v)dxdv. 

Jo "' JJr 2 xS! 



□ 



5.4. An explicit representation formula for M. 
Lemma 7. Let M be the junction defined in US\) . Then 

with the notation 

K = K* ■ • • * K ■ 

n factors 

Proof. Observe that 

r + oo r+oo 

n(t)dt = al e~ at p{t)dt 



(26) 







p(t)e- at dt < 1 , 



where the second equality results from integrating by parts the integral defining k, 
and the final inequality is implied by the fact that p is a C 1 decreasing function. 

By Lemma M S and M > a.e. on R + since f e > a.e. on 

R+ x Z £ x S 1 because f m > a.e. on R 2 xS 1 — see the positivity assumption in 
(fTU| . Applying the Fubini theorem shows that 

M(t)dt= / n(t- s)M(s)dsdt+^ f in (x,v)dxdv n{t)dt 

Jo Jo ™ J v/R 2 xS 1 Jo 

r+oo / r+oo \ p p r+oo 

= M{s) I / K(t - s)dt) ds + 2^// f ln (x,v)dxdv n{t)dt. 

Jo \Js J 7Ta J JrzxS 1 Jo 

In other words 

\\M\\ L , {n+) <\\M\\ Ll0 & +) \\4 L i {n+) + ^ ff f n (x,v)dxdv, 

J JR 2 xS 1 
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so that M € L 1 (R + ) since ||k||l 1 (r+) < 1) an d 

||M|U 1(H+) < _ T , || || r / / f in (x,v)dxdv . 

27rcr(l - ||k|| l i (r+) ) JJ R 2 xS i 

In particular, if 

f m (x,v)dxdv = 



then M = a.e. on R + , so that the representation formula to be established 
obviously holds in this case. 
Otherwise 

// f in {x,v)dxdv > 0; 

define then 

ip(t) := 27TO- f // / in (a:, ujdxcfo') M(t), t > . 

According to Lemma [51 the function ip verifies the integral equation 
(27) ip(t) = (k * (•01n+))(*) + K(t) , a.e. in i > . 

Applying the Fubini theorem as above shows that the linear operator 
A: L 1 (]R + ) 3 f » k* G L 1 (M+) 

satisfies 

/•+0O 

< Mil ll/IU^) with M|| = / «(t)dt < i . 

Therefore (1 — A) is invertible in the class of bounded operators on i 1 (]R + ) with 
inverse 

In particular 

ip = (I — A)' 1 *, = K n 

n>l 

is the unique solution of the integral equation (f2"T)) in L 1 (R + ), which establishes 
the representation formula in the lemma. □ 

5.5. Asymptotic behavior of M in the long time limit. 

5.5.1. The characteristic exponent^. 
Lemma 8. For each a > ; the equation 

o-e-^ + ^p(t)dt = 1 



with unknown £ has a unique real solution This solution £ CT satisfies 

-a < £ a < 0. 

Proof. Consider the Laplace transform of the function n defined above: 

/>oo 

£[«](0 := / ae-^+^p^dt. 
Jo 

As < p < 1, C[k] is of class C 1 on ] — er, +oo[, and 

roc 

£[«](0 = / o-e-^Hpfydt < 
Jo 

as p(t) > for each t > 0. The function is therefore decreasing on ] — cr, +oo[ 
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For each t > 0, 

K,{t)e^^ 0+ as £ -> +00 , 

while 

n{t)e^^ < ae~ at for each t > , 
since < p < 1. By dominated convergence, one concludes that 

£[«](£) -> 0+ as £ -> +00. 

Besides, for each t > 0, 

CTp(i)e- (ff+5)t t op(t) , as £ 4 -a+ . 
By monotone convergence, 

f+oo 

~~ * CT / p(t)dt = +00 , as £ — » — <r + . 
Jo 

(Notice that the equality 

p(t)dt = +00 



follows from the lower bound in (J7J.) 

By the intermediate value theorem, there exists an unique £ CT > —a such that 

£[«](&) = 1- 
Besides £ CT < as £[k] is decreasing and 

C[k](0) = / «(t)d< < / ae- at dt = 1 = £[«](&) , 
Jo Jo 

which concludes the proof. □ 
In particular 

is a decreasing probability density on R + . 

5.5.2. TTie Renewal Equation. It remains to prove statement (3) in Theorem[5] 

First, for each A € R and each locally bounded measurable function / : R i— > R 
supported in R + , denote 

f x (t) := e At /(*) for each t G R. 
Notice that for each such /, g, we have 

e A '(/ * <?)(«) - (fx * ff A )(t) for each * e R. 
Hence, if ip is a solution of the integral equation (f2"T|) , the function ip—£ a satisfies 

(28) v-e„ (*) = * iMJ(t) + , 

which is a renewal integral equation, in the sense of [12] . 

Moreover, as noticed above, K-^ a is a decreasing probability density on R + , so 
that in particular is directly Riemann integrable (see [12] pp. 348-349). Thus, 
applying Theorem 2 on p. 349 in [TJ] shows that 

(29) ^(i)e~ ?CT * -> as t -> +oo. 

By definition of ip, this is precisely the asymptotic behavior of M in Theorem 2 (3). 
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5.6. Two important limiting cases for We conclude our proof of Theorem 
[2] with a discussion of the asymptotic behavior of £ CT (statement (4) of Theorem [5]) 
in the two following regimes: 

(1) the collisionless regime a — > + , and 

(2) the highly collisional regime a — > +oo. 

End of the proof of Theorem^ Denote for the sake of simplicity A CT := a + £ CT . 
Establishing that £ a ~ —a as cr — ► + amounts to proving that X a = o(cr). First, 
notice that, since — a < £ CT , 

< X a < a 

so A (j — > + as cr — * + . Keeping this in mind, we have 

(30) / e- Kt p(t)dt = - 

Jo ° 

by definition of Substituting z = X a t in the integral above, we obtain: 

X f +oa 
< — = / e- z p(z/A CT )dz. 
Jo 

Since A<j — > + as cr — > + and p(t) — > + as £ — » +oo, one has p(z/X a ) — » + as 
cr — » + . Besides < e~ z p(z/\ a ) < e~ z so that, by dominated convergence 

as cr^0+. 

cr 

This establishes the asymptotic behavior of £ CT in the collisionless regime. 

As for the highly collisional regime, we return to the equation (I30[) defining £<j 
(written in terms of A ff ): 

/•+OQ 

1 = cr / e- Kt p(t)dt 
Jo 

/•oo /*oo 

= X a I e^ Kt p{t)dt - £ CT / er Kt p(t)dt 
Jo Jo 

oo 



Jo Jo 
where the last equality follows from integrating by parts the first integral on the 
left hand side. Therefore 



or, after substituting t' = X a t, 
(31) = 



oo 

„-A„t 



p(t)dt 



oo 



oo 

t 



e~ t p{tlX a )dt 
e~*p(i/A CT )df 



Equation ([3"U|) shows that X a — > +oo as cr — » +oo. Passing to the limit in the 
right-hand side of (131|) . we find, by dominated convergence 



e - *p(0)eft 

£<r — * ^joo = p(0) as cr — > +oo . 

e - V(0)eft 
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Indeed p is decreasing and convex, as can be verified for instance on the Boca- 
Zaharescu explicit formula0 (dJ)-® for p, so that 

< -p(t) < -p(0) , for each t > . 

We conclude by observing that the same explicit formulas of Boca-Zaharescu [I] 
imply that 

p(Q) = -2 . 

□ 



6. Final remarks and open problems 

The present work provides a complete description of the homogenization of the 
linear Boltzmann equation for monokinetic particles in the periodic system of holes 
of radius e 2 centered at the vertices of the square lattice el? (Theorem [TJ) In 
particular, we have given an asymptotic equivalent of exponential type of the total 
mass of the particle system in the long time limit (Theorem [SJ) 

Since the discussion in the present paper is restricted to the two dimensional 
setting, it would be useful to extend the results above to the case of higher space 
dimensions, and to lattices other than the square or cubic lattice. Most of the 
arguments considered here can be adapted to these more general cases; however, 
the analogue of the distribution of free path lengths (the function p(t)) is not known 
explicitly so far. See [3] for these more general cases. 

Otherwise, it would also be interesting to investigate other scalings than the 
Boltzmann-Grad type scaling considered here — holes of radius e 2 centered at the 
vertices of a square lattice whose fundamental domain is a square of sise e in the 
case of space dimension 2. Typically, one would like to mix the homogenization 
procedure considered in the present work with the assumption of a highly collisional 
regime d > 1, so that the size of the holes and the distance between neighboring 
holes are scaled in a way that differs from the one considered here. We hope to 
return to this problem in a forthcoming publication. 

Finally, the homogenization result considered in the present paper raises an 
interesting question, of quite general bearing. Usually, homogenization is a limiting 
process leading to a macroscopic description of some material that is known at the 
microscopic scale. In the problem considered here, it has been necessary to use a 
more detailed description of the particle system than that provided by the linear 
Boltzmann equation (problem (5 e ) set in the extended phase space that involves 
the additional variable s.) 

In other words, the formulation of the macroscopic homogenization limit for 
the linear Boltzmann equation considered here involves remnants of an even more 
microscopic description of the system than the linear Boltzmann equation itself — 
namely the extended phase space and the additional variable s. 

We do not know whether this phenomenon (i.e. the need for a more microscopic 
description of a system to arrive at the formulation of a homogenized equation 
for that system) can be observed in homogenization problems other than the one 
considered here — for instance in the case of equations other than those found in 
context of kinetic theory. 



In space dimension higher than 2, one can show that the analogue of p is also nonincrcasing 
and convex, by using a variant of a formula due to L.A. Santalo established in for want of a 
en explicit formula giving the limiting distribution of free path lengths. 
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